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CORRELATION FUNCTIONS FOR RANDOM COMPLEX ZEROES: 
STRONG CLUSTERING AND LOCAL UNIVERSALITY 

F. NAZAROV AND M. SODIN 



Abstract. We prove strong clustering of fc-point correlation functions of zeroes of 
jrt '. Gaussian Entire Functions. In the course of the proof, we also obtain universal local 

bounds for fc-point functions of zeroes of arbitrary nondegenerate Gaussian analytic 

functions. 

rvj , In the second part of the paper, we show that strong clustering yields the asymp- 

^SJ ' totic normality of fluctuations of some linear statistics of zeroes of Gaussian Entire 

Functions, in particular, of the number of zeroes in measurable domains of large area. 

This complements our recent results from the paper "Fluctuations in random complex 

zeroes" . 



1. Introduction 



Consider the Gaussian entire function (G.E.F., for short) F(z) = > Cj—?=j where 

> ; j>o ^ 

C^ , ^j are standard independent Gaussian complex coefficients; that is, the density of the 

probabihty distribution of Q with respect to the Lebesgue measure on the complex plane 
is -e^l'-l . A remarkable feature of the random zero set Zp = F^^{0} is its distribution 

ly-^ . invariance with respect to the isometrics of the complex plane [5, 10, 16]. It's easy to 

^D ' compute the covariance function of F: 

o : ._. 



£{Fiz)Fiw)} 



j_j ■ Hence, after normalization, the covariance equals e™~h\^\''~h\'^\'' = ^^•^{zw)-\\z~wY ^ 

^ which decays very fast when jz — wl grows. This hints at the almost independence 

of random zeroes at large distances. One of convenient ways to formalize the almost 
independence of a zero process at large distances is based on clustering of its fc-point 
functions, cf. [14, § 4.4]. In this paper, we will develop this idea for zeroes of GE.F.'s, 
and then will apply it to prove the asymptotic normality of fluctuations of some linear 
statistics of zeroes of Gaussian Entire Functions. We note that there is a very different 
approach to almost independence of random complex zeroes on large distances that 
proved to be useful in [11, 12, 13]. 

1.1. fc-point functions. The fc-point functions express correlations within fc-point sub- 
sets of the point process. It is customary in statistical mechanics to describe point 
processes by the properties of their fc-point correlation functions. 
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The /c-point function p = pk oi a, random zero process 2 on C is a symmetric function 
p: {Z=(zi,...,Zfc)eC^': Zi^Zj, for i / j} ^ M+ 



defined by the formula 

8{ n #(^ni?,)}= / 



l^i^k 



p{zi, ..., Zk) dA{zi)..AA{zk) 

.xBk 



for any family of mutually disjoint bounded Borel sets Bi, ..., Bk in C. Here A is the 
Lebesgue measure on the complex plane. The /c-point function of zeroes of a Gaussian 
analytic function / exists, provided that for any zi,...,Zk € C with Zi / Zj for i j^ j, 
the random variables /(^i), ■••, f{zk) are linearly independent [5, Corollary 3.4.2]. It is 
not difficult to show (see the beginning of Section 2) that for G.E.F.'s, this condition 
holds for each /c E N. 

1.2. Main results. The first result treats the local behaviour of A:-point functions. It 
appears that for a wide class of non-degenerate Gaussian analytic functions, the A;-point 
functions of their zeroes exhibit universal local repulsion when some of the variables 
zi, ...,Zk approach each other. 

Recall that a Gaussian analytic function (G.a.f., for short) f{z) in a plane domain 
G C C is the sum 

(1) m = Y.Cnfn{z) 

n 

of analytic functions fniz) such that 

{z)\ < CX3 locally uniformly on G, 



Ei/«(^ 



with independent standard complex Gaussian coefficients Cn- 

We postpone until the beginning of Section 2 the technical definition of d-degeneracy, 
which we use in the assumptions of the next theorem. Here, we only mention that 
G.a.f.'s with "deterministic zeroes" (that is, /{zq) = a.s., for some zq G G) are 
1-degenerate. G.a.f.'s such that the random variables /(zi), ..., /(z^) are linearly de- 
pendent for some zi,...,^^ G G, are /c-degenerate, and G.a.f.'s for which the random 
variables /(zi), /'(zi), ..., /(z^), /'(z^) are linearly dependent are 2fc-degenerate. We 
also mention that Gaussian Taylor series (either infinite, or finite) 



fiz) = ^CnCnZ^ 



are d-nondegenerate, provided that cq, ci, ...,Crf_i / 0. In particular, the G.E.F. is 
d-nondegenerate for every positive integer d. 

Theorem 1.1. Let f be a 2k-nondegenerate G.a.f. in a domain G, let pf he a k-point 
function of zeroes of f , and let K d G he a compact set. Then there exists a positive 
constant G = G{k, f, K) such that, for any configuration of pairwise distinct points 
zi,...,Zk G K, 

G~^ JJizi -ZjP ^ pf{zi,...,Zk) ^ GJJIzj - Zjp. 

i<j i<j 
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The next result is a clustering property of zeroes of G.E.F.'s. It says that if the 
variables in C^ can be split into two groups located far from each other, then the 
function p^ almost equals the product of the corresponding factors. This property is 
another manifestation of almost independence of points of the process at large distances. 

For a non-empty subset I = {ii,...,^^} C {1, 2, ..., A;}, we set Zi = {zjj, ..., Zj^}. We 
denote by 

d(Zi,Zj) = inf \zi — Zj\ 

the distance between the configurations Zj and Zj. 

Theorem 1.2. For each k ^ 2, there exist positive constants C^ and A^ such that for 
each configuration Z of size k and each partition of the set of indices {1,2, ...,k} into 
two non-empty subsets I and J with d{Zj,Zj) ^ 2Afc, one has 

(2) 1-5^ /^) <i + ^ with 5 = C,e-|W^-^-')-^^)\ 

p{Zi)p{Zj) 

Combining Theorems 1.1 and 1.2, and taking into account the translation invariance 
of the random zero process Zp, we obtain a uniform estimate for p^ valid in the whole 

Theorem 1.3. For each k ^ 1, there exists a positive constant C^ such that for each 
configuration {zi, ...,Zk), 

Cfc"' n^d^* - ^ii) ^ '^(^i' ■••' ^fc) ^ CkY[i{\zi - zj\) , 

i<j i<j 

where £{t) = min(t^, 1). 

To get Theorem 1.3, we use induction on k. For k = 1 the result is obvious. Now, 
given an integer m ^ 2, suppose that Theorem 1.3 has been proven for fc-point functions 
with k ^ m. To prove it for m + l-point functions, we consider two cases. First, suppose 
that the configuration of points {zi, ...,Zm+i} cannot be split into two groups lying at 
distance 2Am+i or more from each other. Then the result follows from the local bounds 
in Theorem 1.1 and translation invariance. In the other case, the result follows from 
Theorem 1.2 and the inductive assumption. □ 

Note that Theorem 1.3 yields that the fc-point functions pk are uniformly bounded 
on C'^ by constants depending on k. This observation immediately yields the additive 
version of the clustering property: 

Theorem 1.4. For each k ^ 2, there exist positive constants C^ and A^ such that for 
each configuration Z of size k and each partition of the set of indices {1, 2, ..., fc} into 
two non-empty subsets I and J with d{Zj,Zj) ^ 2Afc, one has 

(3) \p{Z) - p{Zi)p{Zj)\ ^ Cfce-^(<^(^^'^-^)-^* 



^kf 



The proofs of Theorems 1.1 and 1.2 start with the classical Kac-Rice-Hammersley 
formula [5, Chapter 3]: 



(4) pf{zi,...,Zk) = / |r/i| ...\r]k\ Vf{r]';zi, ...,Zk)dm{r]i)...dm{r]k), 

Jck 

where Pj( ■ ; zi, ...,Zk) is the density of the joint probability distribution of the random 
variables 

(5) /(zi), /'(zi), ...,/(zfc), /'(zfc). 
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and rj = (0, r/i, ..., 0, r^fc) is a vector in C^ . Since the random variables (5) are complex 
Gaussian, one can rewrite the right-hand side of (4) in a more explicit form 

(6) /3/(2:i,...,Zfc) = 2kA ^j- / l??iP---l^fcpe"2<'^/ '''''>dm(77i)...dm(7?fc), 

vr^'^ det i / Jcfc 

where Fj = rj(zi, ..., Zfc) is the covariance matrix of the random variables (16). We 
consider the linear functionals 

k 

where 






r^{z) 



'' ' "i ^ A 



Z — Zj ' (z — '"'^^ 



and 7 C -fC is a smooth contour that bounds a domain G' <Z K that contains the points 
zi, ..., Zfc. Then we observe that for every vector 6 = (ai, /3i, ..., a^, /3fc) in C^ , we have 

{Tj5,6)=£\Lf\\ 

This observation allows us to estimate the matrix FT , and hence the integral on the 
right-hand side of (6), using some simple tools from the theory of analytic functions 
of one complex variable, and thus avoiding formidable expressions for the Gaussian 
integrals on the right-hand side of (6) that involve quotients of large determinants and 
permanents. 

We note that in [2], Bleher, Shiffman, and Zelditch estimated these large determinants 
and permanents, and proved that if the points Zi are well separated from each other, 
i.e., min \zi — Zj\ ^ r^ > 0, then some estimate similar to (3) holds with a factor C{k, rf) 

on the right-hand side. Unfortunately, in this form their result is difficult to apply. For 
instance, it does not yield boundedness of the A;-point functions on the whole C'^, and 
we could not use it for the proof of the asymptotic normality of some linear statistics, 
see Theorem 1.5 below. 

1.3. Clustering of /c-point functions and asymptotic normality of linear sta- 
tistics. Various ways to derive the asymptotic normality of fluctuations of the number 
of random points in large volumes from clustering of all fc-point functions are known 
is statistical mechanics. Our next result is a variation on this theme. It pertains to 
arbitrary point processes Z on the plane with clustering fc-point functions. 

Let /i be a non-zero bounded measurable function with compact support. For the 
scaled linear statistics n{R\ h) = X^^jg^ ^(■^) ^^ ^^^ point process Z, we set 

n{R- h) = n{R- h) - £{n{R; h)}, a{R; hf = S{n{R; hf], 

In what follows, we consider the functions h with a{R] h) growing as a positive power 
of i?, i.e., satisfying 

(7) a{R; h) ^ R^ 

for some 5 > and for all sufficiently big R. Later, we will see that this holds when h is an 
indicator function of an arbitrary bounded measurable set of positive area (Lemma 1.6 
below) . 
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We call the function ip: (0, oo) — )• (0,oo) fast decreasing if it decreases and for each 
positive m, lim x"^ip{x) = 0. We say that the k-point functions p of a point process Z 

are clustering if there exits a fast decreasing function ip such that for each k ^ 2 and 
for each partition of the set of indices {1, 2, ..., k} into non-empty disjoint subsets / and 
J, one has 

(8) \p{Z) - p{Zj)p{Zj)\ ^ Ck^{ckd{Zi, Zj)) . 

We say that the A;-point functions are bounded, if sup^fc Pk < oo. 

Theorem 1.5. Suppose that the k-point correlation functions pk of a point process 
Z are clustering and bounded. Then for each bounded measurable compactly supported 
function h with cr{R; h) growing as a positive power of R, the normalized linear statistics 
n*{R; h) converge in distribution to the standard Gaussian law as R ^ oo. 

The proof of this theorem is based on estimates for cumulants that follow from clus- 
tering. A similar approach was used by Malyshev in [7], and by Martin and Yalgin in 
[9]. 

We note that there is a counterpart of Theorem 1.5 for arbitrary determinantal point 
processes, namely, a theorem of Soshnikov. In [17], he proved that for arbitrary deter- 
minantal point processes, the fluctuations of linear statistics associated with compactly 
supported bounded positive functions h are asymptotically normal if the variance grows 
at least as a positive power of the expectation. His proof is based on peculiar com- 
binatorial identities for the cumulants of linear statistics that are a special feature of 
determinantal point processes, and is quite different from the one of Theorem 1.5. 

1.4. Lower bounds for the variance of linear statistics. In order to apply Theo- 
rem 1.5, we need to know that the variance of the linear statistics n{R;h) grows as a 
positive power of R. An example of such statistics is the number of points of the process 
Z in a bounded set E of positive measure dilated R times. 

We assume that Z is a translation-invariant point process on the plane with one- and 
two-point correlation functions and that the mean number of points of the process Z 
per unit area equals 1; that is, pi = 1 on M? . We also assume that the 2-point function 
of the point process Z satisfies 

(9) p{zi,Z2) = r{zi - Z2) with r - 1 G L^{M^) 

which is a weak form of clustering. For a set £', we denote by Ie its indicator function. 

Lemma 1.6. Let Z be a translation-invariant random point process on M^ with the 
2-point function satisfying (9). Then there exists a numerical constant c > such that 
for each bounded measurable set E <Z^ of positive area A{E), one has 



a{R, Ie^ ^ cmin ^A{E)R^, ^/ A{E) r\ 



< R< 00. 



Note that we do not impose any smoothness assumption on the boundary of the set 
E. The proof of Lemma 1.6 is based on a simple observation: 

Lemma 1.7. Suppose Z is a translation-invariant random point process on M^ with the 
2-point function satisfying (9). Then there exists a numerical constant C > such that 
for every function h G [L^ n L^)(IR^) and every R> 0, 

(10) aiR,hf^^ f \h{0\'dA{0, 
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where 

h{^) = f h{x)e-^^'''< dA{x) 

is the Fourier transform of h. 

Estimate (10) reduces lower bounds for (t{R, h) to estimates for the tail of the integral 
on the right-hand side of (10). 

1.5. Related results on asymptotic normality of fluctuations of linear statis- 
tics of random complex zeroes. Theorem 1.5 complements our recent results [11] 
on fluctuations in the random complex zeroes Zp obtained by a completely different 
technique. Therein, we showed that if Zp is a zero set of a G.E.F. F{z)^ then the 
following hold: 

(i) If /i is a Co*-function with a > 1, then the fluctuations of n{R, h) are asymptotically 
normal.^ 

(ii) If /i is a Cq function with < a ^ 1, and if for some e > 0, we have u{R, h) ^ i?"""*"^, 
then the fluctuations of n{R, h) are asymptotically normal as well. 

(iii) For each < a < 1, there are C^-functions h with abnormal fluctuations of linear 
statistics n{R, h). 

Acknowledgement. We thank Boris Tsirelson for very useful discussions. 

2. Local universality of /c-point functions. Proof of Theorem 1.1 

Let G C C be a plane domain. We fix a compact set K and a bounded domain Gi 
with a smooth boundary 7 = dGi such that K C Gi and Gi C G. We consider linear 
functionals 

k 

(11) Lg = Y,a,9^"'^Hz,) 

i=i 
with zi,...,Zk S K, acting on functions g analytic in G. Then for every function g 
analytic in G and for every functional L as above, we have 

k 

(12) Lg = l^g{z)r^{z)dz with r^(^) = ^ E (J^fj^ ' 

Note that r^ is a rational function vanishing at infinity. We put 

rank(L) = deg(r'^) . 

Then we say that a G.a.f. f{z) in G is d-degenerate, if there exists a functional L of 
rank at most d such that 

(13) Lf = almost surely. 

Otherwise, we say that / is d-nondegenerate. For instance, G.a.f. 's with "deterministic 
zeroes" (that is, f{zo) = a.s., for some zq G G) are 1-degenerate. G.a.f. 's such that the 
random variables /(zi), ..., /(z^) are linearly dependent are /c-degenerate, and G.a.f.'s 
for which the random variables /(zi), /'(-?i), •••, fiz^), f'{zk) are linearly dependent are 
2fe-degenerate. 



This result was preceded by yet difTerent proofs in [16, 18] of the asymptotic normality in the case 
of Co-functions h. 
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Observe that the G.a.f. (1) satisfies condition (13) if and only if L/„ = for every n. 
So Gaussian Taylor series (infinite, or finite) 



fiz) = ^CnCnZ^ 



are d-nondegenerate provided that cq, ci, ..., Cd„i ^ 0. Indeed, suppose that (13) holds. 
Then 

Lz"= f z''r^{z)dz = 0, for n = 0,l,...,d-l. 

Choose R so large that G C {\z\ < R}. Then 

z"r^(z) dz= z"r^(z) dz = 0, for n = 0, 1, ...,d - 1 , 

{\z\=R} A 

whence, 

r-^(z) = 0(z-('^+i)), z^oo. 

That is, the degree of r^ is not less than d+1, which means that f{z) is d-nondegenerate, 
as we claimed. In particular, we see that the G.E.F. is d-nondegenerate for every non- 
negative integer d. 

Claim 2.1. Let f be a d-nondegenerate G.a.f., and let K C G be a compact set. Then 
there exists positive constants c{d, f, K) and G{f, K) such that for each functional L of 
rank at most d with zi, ..., z/;. £ K, we have 

(14) cmax|r^p ^f|L/p ^Cmax|r^|2. 

7 7 

Proof of the upper bound in (14); 

2 



£\Lf\' = £ 



/r^dz 

7 



2 p / \ 2 



^ Aength(7) • max Ir^l") • j £\f\^ \dz\ = C(/, 7) T: 



max Ir I I . □ 



7 V ^ 

The proof of the lower bound uses a simple compactness argument. Denote by Rd{K) 
the set of rational functions of degree at most d vanishing at 00 and having all poles in 
K. 

Claim 2.2. Each sequence of rational functions {rm} C Rd{K) with max^ {r^l ^ 1 
has a subsequence that converges uniformly on j to a function r G Rd[K). 

Proof of Claim 2.2: We have r-m = 2^, where Pm and qm are polynomials of degree 

^d—1 and ^ d correspondingly, and ^^(-2) = Wi=i{z — Zi{m)) with zi(m), ..., Zrf(m) GiC 
for each m. Choosing a subsequence, we may assume that q^ converge uniformly on 
7 to a polynomial q of the same form. Then max^ \Pm\ ^ C for each m. Since the 
degrees of the polynomials pm are uniformly bounded, we can choose a subsequence 
that converges uniformly on 7 to a polynomial of degree ^ d — 1. □ 

Proof of the lower bound in (14).- Suppose that there exists a sequence of functionals 
Lm, of the form (11) of rank d or less, with max^ |j-im| — \ ^^^ with £\Ljnf\'^ — ?• 0. 
By (12) and Claim 2.2, we can choose a subsequence of this sequence such that r^™ 
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converge uniformly on 7 to some nonzero r E Rd{K) (to simplify our notation, we omit 
subindices) . Then r = r for some functional L of the same form (11), and 



hm £\Lmf - Lf\^ ^ (length(7))2 • [ £\f\^\dz\ ■ lim 



maxlr^™ - r^P = 0. 



¥00 7 

Hence, Lf = almost surely, which contradicts the d-nondegeneracy of /. □ 

To estimate the A;-point function, we use a classical formula that goes back to Kac, 
Rice, and Hammersley, see [2] and [5, Chapter 3]: 

(15) pf{z,,...,Zk)= ,,.\^ / |r?i|2...|r/fc|V5<r7V,V)d^(^^)...dA(7?fc), 

7r^'= det Ff J^k 

where Fj = rj(zi, ..., z^) is the covariance matrix of the random variables 

(16) /(zi),/'(zi),...,/(zfc),/'(z,), 

and r]' = (0,r/i, ..., 0,7^^) is a vector in C . Here, we assume that the random vari- 
ables (16) are linearly independent; later, we will impose on the G.a.f. f{z) a somewhat 
stronger restriction of non-degeneracy. The Kac-Rice-Hammersley formula (15) allows 
us to reduce estimates for the A;-point function pf to estimates for the covariance matrix 

We start with a simple observation. Consider "special" linear functionals of rank 2k: 

k 

(17) L/ = j;[a,/(z,) + /3,/'(^,)]- 

Claim 2.3. For every vector 5 = (ai,/3i, ...,ak, Pk)"^ in C?^ , we have (T/iJ, 5) = E\Lf\^ . 
Proof: by straightforward inspection. D 

Introduce the Gaussian polynomial 



f2k-l{z) = ^Cn 



2k-l 
n=0 



Claim 2.4. Let f be a 2k-nondegenerate G.a.f. and let K G G be a compact set. Then 
there exists a positive constant C{k, f, K) such that for each configuration zi, ..., Zk G K, 

(18) C-^Tf^^_^{zi,...,Zk) ^Tf{zi,...,zk) ^ CTf^^_^{zx,...,Zk), 
where the inequalities are understood in the operator sense. 

Proof: this is a straightforward consequence of Claims 2.3 and 2.1. D 

Now, using the Kac-Rice-Hammersley formula (15), we readily get equivalence of the 
/c-point functions pj and Pf2^_^'- 

Claim 2.5. Let f be a 2k-nondegenerate G.a.f. and let K G G be a compact set. Then 
there exists a positive constant C{k, /, K) such that for each configuration zi, ..., Zk G K , 

(19) C~^pf2^_^{zi,...,Zk) ^pf{zi,...,Zk) ^ Cpf2^_^{zi,...,Zk). 
Proof: First, note that Claim 2.4 yields 

(20) C-'" det F;,,„, ^ det F; ^ C" det F^..^, . 
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Then, plugging estimates (18) and (20) into (15), we get 

C"' f ,.,2 ,..,2 -|c-(r- VV) 



Pfiz,,...,Zk)^^j—— / \m\'-\Vk\'e -'-- ^^^2^-i'""d^(77i)...dA(%) 

TT^'^detT/j^^^ Jck 



C 



2k 



■C'^ / \vi\'...\m\'e--^^f^^-^'''''^dA{v,)...dA{r„) 



7r2'= det r/2,_^ Jck 

Similarly, we get the lower bound in (19). □ 

We see that it suffices to estimate the /c-point function only for one special Gauss- 
ian polynomial /2fc-i of degree 2k — 1. First we consider the probability density 
Pf2^._i{zi, ..., Z2k-i) of the joint distribution of all zeroes of f2k-i- This means that 
Pf2k-i ^^ ^ symmetric function such that, for any symmetric function S* of 2A; — 1 com- 
plex variables, the expected value of S{zi, ..., Z2k-i) equals 

/ S{zi, ..., Z2k-l)p2k-l{zi, ..., Z2k~-l) d^(2;i)...d^(z2fc-i) • 

The density p/j^.i can be computed using a classical formula for the Jacobian of the 
transformation of zeroes of the polynomial into its Taylor coefficients [3, 4]: 

p/2fc-i(^i'-'^2fc-i) =Cfc Yi \^i-^i?-{ 5Z l^Jf)~^^ 

where Oj 's are the coefficients of the polynomial 

JJ (z- Zi)= ^ OjZ^ (here, (J2k-\ = 1)- 

l^j^2fe-l 0^js£2fc-l 

That is, 

P/2fc-l(2;i'---'^2fc-l) = W \Zi- Zjl"^ ■ H{zi,...,Z2k-l), 
l^i<j^k 

where H is a non-negative continuous function. 

Claim 2.6. There exists a positive constant C{k, K) such that for all zi, ...,Zk€K , and 
all Zk+i, ..., Z2k-i S C, we have 

2k~l 

H{z,,...,Z2k^i)^C{k,K) II {l + \zi\)-^. 

i=k+l 

Proof of Claim 2. 6: Let 

2fc-l 

P{z)= J{{z-Zi). 

i=l 

For each z G T, we have 

2fc-l 2fc-l 

i=o i=o 
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An elementary geometry shows that there exists a point z £ T so that \z — Zi\ ^ Ck for 
z = 1, ..., 2k — 1, whence, \z — Zi\ ^ Cfc(l + \zi\) for i = 1, ...,2k — 1, and 

2k-l 2k-l 



•J=A:+1 



Therefore, 



and 



2fe-l 2fc-l 

j=0 j=fc+l 



2A:-1 

4fc 



At last, using that \zi — Zj\ ^ (1 + |zi|)(l + \zj\), we get 

2A:-1 2fc-l 

/7(zi,...,Z2fc-i)^C(A:,/0 n (l + l^^l)'^''"'^ 11(1 + 1^^1)""" 

i=fc+l i=k+l 

2fe-l 

=c(fe,Jo n (i+i^ii)"' 

i=k+l 

completing the proof of the claim. D 

Now, we readily finish the proof of the local universality theorem. The fe-point func- 
tion pf^^_^{zi,....,Zk) can be obtained from the probability density function p/2fe_i by 
integrating out the extra variables 



/'/2fc-i(^l'"-'^fc) = '^k I P/2fc_i(^l,---,^2fc-l)dA(^fe+l)---d^(22fc-l) 

= Ck Y{ l^*-^il^"/ H{zi,...,Z2k-i)dA{zk+i)...dA{z2k~i). 

Due to Claim 2.6, the integral on the right-hand side converges uniformly with respect 
to the variables zi,...,Zk E K. Hence, it is a continuous function in zi, ..., Zk. Given 
zi, ..., Zk, the function H (as a function of variables Zk+i, ..., Z2k-i) vanishes only on 
a set of zero Lebesgue measure in C , hence, the integral is a positive function of zi, 
..., Zk. This completes the proof of Theorem 1.1. D 

3. Clustering of /c-point functions. Proof of Theorem 1.2 

The proof again employs the linear functionals L introduced in (17) and Kac-Rice- 
Hammersley's formula (15). By pT we denote the circumference {\z\ = /?}. 
Let i^ be a G.E.F. . For w £ C, denote by Tyj the projective shift 

T^F{z) = F{w + z)e-^^e-5l'"l' . 

Comparing the covariances £ \ F{zi)F{z2) > and £ I TwF{zi)T^F{z2) >, one readily checks 
that for every w £ C, T^F is also a G.E.F. . Our first aim is to estimate the covariance 
£ \Li [Tyj^F) L2 {Tyj^F) >, where Li, L2 are two linear functionals of the kind (17), and 

Wl,W2 eC, Wi y^ W2. 
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Claim 3.1. Let k ^ 1, p > 0, \wi — W2\ ^ 4p, and let the points Zj participating in the 
definition of the functional L satisfy \zj\ ^ p (j = 1, ..., k). Then 

(21) |£:{Li(T^,F)L2(T^2^)} 

Note that, since T^F is equidistributed with F for aW w £ C, we can also rewrite the 
sum on the right-hand side as S\LiF\'^ + S\L2F\'^. 

Proof of Claim 3.1: Using the representation (12) with 7 = 2pT, we can write 
£ [Li{T^,F)L2{T^,F)^ 

If r^^{zi)rL^{z2)£{T^,F{zi)T.^,F{z2)\dziAz2 

J J2pJx2pJ ^ ' 



2pTx2pTI 



(maxlr ^1 -|-max|r ^1 ) • sup 

^ 2pT 2pT ^ 2i,^2e2pT 



£ [T^,F{zr)T^,F{z2)] 



Since for any «;i, tt;2, Ai, A2 G C, we have 



£[t^,F{\i)T^,F{\2)} 



62 



5|Ai|2 + i|A2|2-i|(t«i+Ai)-(«;2+A2)|' 



the last supremum does not exceed e^P e zd^'i "'2 1 4p) ^ provided that |wi — W2I ^ 4p. 
Taking into account that |r^p ^ C(A;, p)(5|L.jFp everywhere on 2pT (Claim 2.1), we 
get (21). 

In what follows, we also need the following simple geometric claim. 



D 



Claim 3.2. Suppose that zi,...,Zk G C, and d: (0,-|-oo) — )■ (0, +00) is any increasing 
function. Then there exists p ^ I and a covering of points Zj (j = l,...,k) by disks 
D{wn,p) (n = 1,...,N) such that N ^ k and \wn' — Wn"\ ^ d{p) for all n' ^ n" . 
Moreover, p is bounded from above by some constant depending on d and k only. 

Proof of Claim 3.2: Take pi = 1 and consider the covering by the disks D{wn, Pi) with 
N = k, Wn = Zn for all n. If all pairwise distances \wn' — Wn"\ ^ d{pi), we are done. 
If not, replace two centers Wn' and Wn" satisfying \wn' — Wn"\ < d{pi) by one center 
\ {wn' + Wn" ) and increase all radii to P2 = Pi + ^d{pi ) . If in the resulting configuration 
all pairwise distances between centers are at least d{p2), we are done. Otherwise, again, 
replace two close centers by one and increase the radii to ps = P2 + \d{p2), and so on. 
The process will stop after at most k — 1 steps, and the final radius will be bounded 
by the A;-th term in the recursive sequence given by pi = 1, and pj+i = pj + \d{pj). □ 

Claim 3.3. There exists a constant A = A(/c) with the following property. Let Zi, ..., z^ £ 
C. Suppose that the set of indices {1,2, ...,k} is partitioned into two non-empty subsets 
I and J so that 

d = d{Zi,Zj) = mi{\zi-Zj\: i£L,j eJ}^ A. 

Then for any two linear functionals 

L'f = 5][a./(z.) + ft/'(z.)], L'f = J][a,/(z,) + /3,/'(z,)], 
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(22) 



£[{L'F){LJF)}\^C{k)e--^^''-'^^'{£\L'Ff+£\L^F\^} . 



Proof of Claim 3.3: We put 



d{p) = Ap+ y/2\og{lQkC{p,k)) 

where the constant C{p, k) is the same as in (21), and apply the construction of Claim 3.2 
to the points Zj. We get a number p ^ p{k) with p{k) ^ 1 independent of zi, ..., z^, and 
a sequence of disjoint disks D{wn,p) {n = 1, ■■■,N, N ^ k) such that each point Zj lies 
in at least one of these disks. Then we take A = 6p{k), and notice that no disk can 
contain two points Zi, Zj with i £ I, and j G J. Thus, the set of indices {1,2, ..., A^} 
gets partitioned into two non-empty subsets I and J such that the disks corresponding 
to indices from I cover the points Zi with indices i G I, and similarly for J and J7. 
Now, we put 

Lnf= Yl h/(^0 + ft/'(^i)] 

i: ZiGD{w„,p) 

and note that Lnf = {Ln o T^^^ C^Wuf) = Ln {Tw„f), where L„ is some linear functional 
of the kind (17) in whose definition the points Zi — Wn with Zi G D{wn,p) participate ^. 
Observe that 

£\L'F\^ = Y,£\LnT^y^ 

neX 






{(L„,r^„,F)(L„»r^,,,F)} 



=a{n' ,n") 



Y,s\LnF + Yl <n\n"). 



By Claim 3.1, for n' ^ n", we have 

\a{n',n")\ ^ C(/>, A:)e-5(^(^)-4p)' 



£ 



Ln'F 



neX 



+ £ 



n',n"eX,n'=/=n" 



Ln"F 



< 



mk 



£ 



Ln'F 



+ £ 



Ln"F 



by our choice of d{p) , so 



I I 9 ^—^ 



n€X 



A similar estimate holds for L , whence. 



£ \L^F\ + £ iL-^Fl ^ -Y£\LnF 

I I ' I I ^ 9 Z^ " 



On the other hand. 



£{{L^F){yF)]\^ Y Hn',n")\. 



n'eX,n"eJ 



^More explicitly, if Lf = E, [a^fiz^) + |3^f'{z^)] , then 

i : 2 j£Z){tD,j ,p) 
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I Wn' — Wn" I ^ \zi — Zj\ — 2p ^ d — 2p . 



where Zi is any point in D{'Wni,p), and Zj is any point in D{wn",p). Using Claim 3.1 
again, we get 



|a(n',n")| ^ C(y9, A:)e-5('^-6p)' 



£ 



Lr,lF 



+ £ 



L/ri"F 



and 



J2 \<^in\n")\ ^ C{p,k)e~'2(d-^P^' J^^l^nF 

n'£X,n"£j n 

Recalling that p ^ p{k), we get the claim. 



n 



Now, we are ready to prove Theorem 1.2, that is, the multiplicative clustering. To 
simplify our notation, we denote by 

e = e(d,A;) = C(A;)e-5('^-^)' 

the small factor on the right-hand side of (22). Recall that d is a shortcut for d{Zj, Zj) 
and that d ^ 2 A. In what follows, we assume that A is so large that e < g- We will 
show that 



(23) 



l-£ 

1 + e 



-I 2k 



p{Zi)p{Zj) ^ p{Z) ^ 



l + e 
1-e 



-I 2k 



p{Zi)p{Zj) , 



which yields Theorem 1.2. 

To prove (23), we repeat the argument already used in the proof of Theorem 1.1, 
and interpret inequality (22) from the previous claim as a two-sided estimate for the 
covariance matrix T = Tp{zi, ...,Zk) of the family of 2k Gaussian random variables 
{F{zi),F'{zi)) (1 ^ i ^ k). Consider also the covariance matrices Tj and Tj of the 
subfamilies with i £ I and i £ J respectively. Then estimate (22) can be rewritten in 
the following form 



(24) 
where 



(i-e)rj,j^r^(i + e)r,,j 



r/„ 



Ti 

r. 



is the covariance matrix corresponding to two independent copies Fj and Fj of the 
G.E.F. F. Inequalities (24) yield 

(1 - ef'' det Tij ^ det F ^ (1 + ef'' det Tij . 

We again use the Kac-Rice-Hammersley formula (15): 

1 



p{zi,...,Zk) 



tt"^^ det r 



m\' -IVkl' e-"^^"^-"^''^'^ dA{m) ...dAir^k) 



where r]' = (0, r]i, ..., 0, rjk) is a vector in C^ . Note now that 

(l + e)-2'=(detr7)-^(detrj)"i ^ (detT)"^ ^ (1 - e)"^^' (detT/)"^ (det Tj)"^ , 
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and that ir = vr^* • vr^-', where i = \I\ and j = \J\. To compare the integrals, observe 
that 

/ |r?ip ... |%p e-i<r-^VV) dA(r/i) ... d^(%) 

= {l-ef^ [ |r?i|2...|%|2e-5<r7>''''')dA(r/i)...d^(%) 

= (!-£)''= / |r/i|2...|r/,|2e-^<r"''^'''^>d^(r/i)...dA(77,) 

|?7i|2 ... |r?jf e-^<rjS:„Vj> d^(77i) ... dA{r]j) , 

where ??' = (0,r/i, ..., 0,%)^, ??j = (0,r/i, ...,0,77^)^, and r/^ = (0,r/i, ..., 0,7/^)^. Similarly, 
we get an upper bound of the same kind with the factor (1 + e) . This proves (23) and 
hence proves Theorem 1.2. D 

4. Central limit theorem via clustering and cumulants. 
Proof of Theorem 1.5 

The proof of Theorem 1.5 is based on the estimates of cumulants Sk = Sk{h) of the 
linear statistics n{h) = J2a€Z ^i^) defined by the formal expansion 

iog£:{eM'^)} = E^^'=- 

We will show that for each bounded measurable function h with compact support spt(/i) 
and for each k ^ 2, 

\sk{h)\^CkA{spt{h))\\h\t. 

The key point of this estimate is that applying it to the scaled statistics n{R; h) we 
get the factor R^ on the right-hand side. Together with the growth of the standard 
deviation (t{R; h) as a positive power of R, this will yield the asymptotic normality of 
n{R] h) when R — )• 00. 

First, we introduce some notation and recall well-known facts about the moments and 
cumulants. We denote by I\.{k,j) the collection of all unordered partitions of the set 
{1, 2, ..., k} into j disjoint nonempty blocks, and by n(/c) the collection of all unordered 
partitions of the set {1, 2, ..., k} into disjoint nonempty blocks. For vr S n(A;), we denote 
the blocks by {vri, ...,t^j} with an arbitrarily chosen enumeration, and denote the lengths 
of the blocks by pt = |7rt|, 1 ^ t ^ j. The number of blocks in the partition vr will be 
denoted by |7r|. 

The moments rii^ and the cumulants s^ are related to each other by classical formu- 
las [15, Chapter II, § 12]: 

k 

(25) ,, = J^(_i/-i(^_l)! Y, 

i=i TT£U{k,e) 

The inverse formula has the form 

(26) mk= Yl 

7rgn(fc) 
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Next, we recall the definition of so called truncated (or Ursell) /c-point functions p]^. 
The truncated /c-point functions play the same role for the cumulants as the usual 
A:-point functions for the moments. Their definition is suggested by (25). They are 
symmetric functions 

pI: {Z = {zi,...,Zk)£C'': Zi^Zj, for i/j} ^M 

defined by the formula 

k 

(27) pT(^z) = Y,{-^Y"Hi-iy- E Pn{Zn,)-.Pp,{Z^,). 

i=i 7ren(fc/) 

In particular, pj = pi, 

(28) p^{zi,Z2) = P2{Z1,Z2) - pi{zi)pi{z2) , 

then 

pI{zi,Z2,Z3) = P3{ZI,Z2,Z3) 

- {P2{zi,Z2)piiz-i) + P2{zi,Z-s)pi{z2) +^2(^2,^3)^1(^1)) 

+ 2pi{zi)pi{z2)pi{z3), 

and so on .... The inversions to (27) look the same as in (26): 

(29) Pk{Z)= Y. PpAZn,).:pl{Z.J. 

7ren(fc) 

See e.g. [8, Appendix A. 7] for the proof of the equivalence of the formulas (27) and (29) 
based on generating functions; another short and direct proof of this equivalence can be 
found in [1, § 2]. 

Asymptotic factorization of /c-point correlation functions expressed in (8) yields fast 
asymptotic decay of truncated A;-point functions when the diameter of the configuration 
gets large. 

Claim 4.1. Suppose that the k-point functions pk of a point process Z are clustering, 
i.e., satisfy (8) with a fast decreasing function ip, and are bounded. Then, 

(30) |p^(Z)| ^Cfe^(cfcdiam(Z)), 

where ip = m.m{l,ip) and diam(Z) = max \zi — zA is the diameter of the configura- 
tion Z . 

Proof of Claim 4-1'- We show that for each /c ^ 2 and for each partition of the set of 
indices {1, 2, ..., k} = I U J, 

(31) \pl{Z)\^Ck^{ckd{Zi,Zj)). 

Since for each configuration Z, there exists a partition {1, 2, ..., k} = lUJ with d{Zj, Zj) ^ 
Cfcdiam(Z), estimate (31) yields (30). 

To prove (31), we use induction on k. For /c = 2, estimate (31) follows from equa- 
tion (28), clustering of p2, and boundedness of pi and p2. Now, suppose that fc ^ 3. 
We fix a partition {1, 2, ..., k} = I \J J with i = \I\ and j = \J\. We say that a partition 
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vr E n(fc) refines the partition /U J if each set tTs is a subset of either /, or J. Otherwise, 
we say that the partition vr mixes I U J. By the inversion formula (29), we have 

PkiZ) = pk{Z)- Yl Ppii^^i)--- Ppt i^^t), 

t=\Tc\^2 

and 

Pi{Zi)pj{Zj) = ^ pI^ {Z^^ )--Ppt {Znt ) • 

IT refines lUj, t=\n\^2 

Therefore, 

pUz) = pk{z) - pi{Zj)p,{Zj) + Y. pI,{z^,)...pI{z^,). 

IT mixes lUj, t=|7r|^2 

By the induction assumption, each term in the sum on the RHS contains at least one 
factor which is bounded by the right-hand side of (31), while the other factors are 
bounded by constants. □ 

Claim 4.2. Suppose that the truncated k-point functions pj. of a point process Z sat- 
isfy (30) with a fast decreasing function if. Then, for each k ^ 2, 

(32) sup/ \pl{Z)\dA{zi)...dA{zk^i)<oo. 

Proof of Claim 4-2: We fix a point z^ G C, and set Z = [Z',Zk), where Z' G C'^~^. 
Then we split C'^~^ into disjoint sets Gf. 

Go = {Z': diam(Z) s^ 1} , and d = iz': 2^"^ < diam(Z) ^ 2^\ for i^l. 

Applying estimate (30), we get \pj^{Z) | ^ Ck^ (cfc2^~^) provided that Z' G Gi with £^1. 
Then the integral of the absolute value of pjl over Gg, does not exceed 

C,22^(^-i)^(c,2^"i). 

Hence, the integral on the left-hand side of (32) does not exceed 

C, + C,J]22^(^-i)^(c,2^-i)<C, 

(in the last estimate we use that the function ip decreases faster than any power). □ 

The computation of the /c-th cumulant Sk{h) of the linear statistics n[h) requires 
knowledge of j-point truncated functions with all j ^ k. For bounded compactly sup- 
ported measurable functions hi, ..., hk on C, we put 

{hi,...,hk)'^ = / hi{zi) ...hkizk)p'^izi,...,Zk)dA{zi)...dA{zk) . 



/c'= 
Given a set P = {pi, ■.■,Pk} of k positive integers, we define 

s{h,P) = {hP\...,hP^f. 

Given a partition 7, we denote by P^ = {|7i|, ..., |7j|} the set whose elements are the 
lengths of the blocks of 7. 

Claim 4.3. We have 

(33) Sk{h)= Y <h,P^). 

7en(fc) 
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Probably, these relations are well known to those who worked with fc-point functions 
of point processes. For this reason, we relegate the proof of Claim 4.3 to the appendix. 
Now, we easily estimate the cumulants of the linear statistics n(/i): 

Claim 4.4. For each hounded compactly supported function h, we have 

\sk{h)\ <.Ck\\h\t,A{spt{h)). 

Proof of Claim 4 ■4'- By equation (33), it suffices to estimate from above |(/i^S •••; hP^)'^\ 
where pi, ...iPj are positive integers such that ^jPj = k. 
Let G = spt(/i). Then 

\{hP\...,hP^f\ ^ f \h{z^)r ...\h{z,)\P^ \p'^{z^,...,Zj)\dA{zi)...dA{zj) 

JO 

^\\h\L I M^l)-M^j)\p'^i^l,-,^,)\dA{z^)...dA{z,) 

Jo 
^ \\h\t I tG{z,)dA{z,) [ |p^(zi, ...zj)\dA{zi)...dA{z,^,) 

? Ck\\h\t I MzMMzj) = CkMt A{G) , 

JC 
proving the claim. D 

Proof of Theorem 1.5: Denoting by Sk{R\ h) the cumulants of n{R; h), and by s^(i?, h) 
the cumulants of the normalized random variables n*{R;h), we have s^ = 0, ^2 = 1, 
and for k ^ 3, 

*(n h\ ^k{R]h) 

By the estimate from the previous claim, we have 

\sk{R;h)\ i^C{k,h)R^ . 
Recalling that a{R] h) grows as a power of R, we see that for large enough /s's, 

lim sl{R\h) =0. 

By a version of the classical theorem of Marcinkiewicz (see [6, 17]), this suffices to 
conclude that the random variables n*{R;h) converge in distribution to the Gaussian 
law when ii — )• oo. This finishes off the proof of Theorem 1.5. □ 

5. Lower bound for the variance. Proof of Lemmas 1.6 and 1.7 
First, we prove Lemma 1.7. 

Proof of Lemma 1.7: It suffices to prove the lemma for R = 1. The general case readily 
follows by scaling 

x^Rx, h^h{R-^-), C^R-^C, h^R^h{R-). 
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Recall that p{zi,Z2) = r{zi — Z2), and let k,{z) = r{z) — 1. By the assumptions of the 
lemma, this is an L^(M^)-function. We have 



(7(1; hy= I h{zi)h{z2) [p{zi,Z2) + S{ZI - Z2) - 1] dA{zi)dA{z2) 

h{zi)h{z2) [k{zi - Z2) + 5{zi - Z2)] dA{zi)dA{z2) 



(in the last line we used Parseval's formula). By the Riemann-Lebesgue lemma, k(,^) — )■ 
when 1^1 — )• 00. Thus, choosing a sufficiently big constant C so that |k(^)| ^ ^ for 
lei ^ C, we get 






proving the lemma. □ 

Remark 5.1 (cf. Martin- Yalgin [9]). The computation above gives us 

= [1 + k(0) + o(1)] \\h\\l,^^,^R\ 

as i? — )• 00. We conclude that translation-invariant point processes with r — 1 G L^(M?) 
can be divided into two groups: the processes with the Poissonian behaviour of the 
variance of linear statistics, when k(0) 7^ —1, i.e.. 



/ [r(x) - 1] dA(x) / -1 

JR2 



and the processes with non-Poissonian behaviour of the variance of linear statistics, 
when k(0) = —1, i.e., 

(34) / [r{x) - 1] dA{x) = -1 . 

Jr2 

Sometimes, the latter processes are called "superhomogeneous processes". Lemma 1.6 
proven below is trivial for the processes from the first group. On the other hand, 
though condition (34) looks like a degeneration, many interesting translation-invariant 
processes including our main hero, the zero point process Zp, and particle processes 
based on Coulomb interaction that occur in statistical mechanics belong to the second 
group. 

Remark 5.2. For the random zero process Zp, the computation started above can be 
continued. This way, we arrive at an explicit formula for the variance [10] which yields 
that 

a{R;hf^R-^ f \hiO\'\^\UA{0 + R' f |MO|'d^(0- 

Now, we prove the remaining Lemma 1.6. 

Proof of Lemma 1.6: Fix a multiplier m € C^{^) with the following properties: 
^ m ^ 1 everywhere, 'm{^) = for |^| ^ 2, and m{^) = 1 for |^| ^ 1, and set 
m„(^) = m{R~^^). We use the function m/j to cut high-frequency oscillations in the 
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spectrum of the indicator function !„. Denote by ip^ „ the inverse Fourier transform of 

the product 1^ "'"^r- Since VmZ{x) = R^{Vfh){Rx), we have ||Vm^||^i = i?||Vm||2,i, 
whence 

(35) \^^e,r\^R- 

By Lemma 1.7 we have 

a{R,tj^f>R^ f fy^^R^ f \r^\'\l-m \^ = R^ f \1^ - ip 
J\ii\^R' Jv " Jr^ 

with R' = CR, where C is a constant on the right-hand side of estimate (10) in the 
statement of Lemma 1.7. We estimate the integral on right-hand side twice. The first 

bound works when A [{y^p r' ^ |}) ^ |^(-^)- 

W-^E,R'\"^ i J\-^E,R'\" 



E,R 



1 1 5 



1 f ^-^^l 1 f 

4 A <n^ ^ Rh <u'^^^.«'' 

1 /-l/^ / N 

/ Length f {ip , = t}\ dt (coarea formula) 

1 /-i/a I -, T 

-— / \ATiA{{ip, ^ t} j dt (isoperimetric inequality) 



R 



^^■^(Kh'HO^^v^M. 



On the other hand, if A ( {if r/ ^ ^} ) ^ \-^{^)i ^^ have 



/ 



\^E-^E,R'\'^ I ^\^E-^E,R'\" 






In both cases, 



R^ [ fy^> mm \A{E)R^y^A(E)R} 



proving the lemma. □ 

Appendix: Proof of Claim 4.3 

First, we prove relations analogous to (33) that express the A:-th moment mk{h) of 
the linear statistics n{h) in terms of j-point functions pj with j ^ k. For bounded 
compactly supported functions hi, ..., /ijt on C, we put 

{hi,...,hk) = / hi{zi) ...hkizk) pizl,...,Zk)dA{zl)...dA{zk)■ 
Given a set P = {pi, ■.■,Pk} of k positive integers, we define 

m{h,P) = {hP\...,hF^). 
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As above, given a partition 7, we denote by P-y = {|7i|, ..., |7j|} the set whose elements 
are the lengths of the blocks of 7. Then 

(36) mk{h) = y^ m{h,P^). 

7Gn(fc) 

To prove this relation, we denote by (oi, ...,afc) an ordered finite sequence of length 
k of points in Z, possibly with repetitions. Then 



mk{h)=£<. ^ h{ai)...h{ak] 



_(ai,...,afc) 

where the sum on the right-hand side is taken over all finite sequences (ai, ..., afc)- We 
fix a partition 7 G n(fc) and say that the sequence (ai, ..., a^) is subordinated to 7 if the 
following condition holds: a^ = at for s ^ t if and only if the indices s and t belong to 
the same block of the partition 7. In this case, we'll write (ai, ...,a/c) < 7. Clearly, each 
finite sequence of length k is subordinated to one and only one partition 7 G n(/c). 

Given a partition 7 € Il{k,j), we fix an arbitrary enumeration 71, ...,7^- of its blocks. 
Then we have 



^ h{ai)...h{ak) = 

(ai,...,afe)<l7 {bi,...,bj) 

bs^bt for sj^t 



J2 hM{b^)...h\^^\{b,), 



and then 



(ai,...,afc) 



i=i 7en(fcj) 



5] /il^^l(6i).../^l^^l(6,0 



(fei,...,6,) 
bsT^fct for s^i 



> . 



Since 



£\ Y. h,ibi)...h,ib,) 

(bi,...,bj) 
^bs^bt for sj^t 



{hi,...,hj) 



D 



(see for instance, [5, Section 1.2]), we complete the proof of relations (36). 

Now, we turn to the proof of relations (33). Given a partition 7 G Il{k,j) we fix an 
arbitrary enumeration 71, ..., jj of its blocks. Then according to the definition (27) of 
truncated functions p , we have 



(37) s{h,P^) = Y,{ 



-iY-\e 



I) 



E 

neU{j,t 



C3 



h\^'\{zi)...h\^^\{z,)p{Z^,)...p{ZnJdA{zi)...dAiz,) 



where the partition vr splits the set of blocks {71,..., 7^} into i 'super-blocks' H*, 1 ^ 
t ^ i, which we also enumerate arbitrarily. By pt, 1 ^ t ^ i, we denote the number of 
blocks of 7 that are included into the super-block 11*. By P* = P* = | |7qi |, ..., \jap^ 1} 
we denote the set of the lengths of the blocks 7q,^, ..., 70,^^ that are included into the 
super-block H*, this is the set oi pt positive integers. For instance, let A; = 8, j = 3, 
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1 = 2, let 71 = {1, 2, 3}, 72 = {4, 5, 6}, 73 = {7, 8}, and let vri = {1, 2}, ■K2 = {3}. Then 
ni = {71,72}, n2 = {73}, pi = 2, p2 = 1, and pi = {3,3}, P^ = {2}. 

We factor the integrals on the right-hand side of (37) 

I 
/ h\'^^\{zi)...h\^^\{z^)p{Z^,)...p{Z^,)dA{zi)..AA{z,) = \\m{h,P'^^;) . 

Ja fJi 

Then the right-hand side of (33) equals 

(38) E E E(-i)'"^(^-i)' E n-(^4J- 

i=i 7en(fcj) £=i ireTiu/) t=i 

We say that a partition 7 G Il(k,j) refines partition a G Il{k,i), i ^ j, if each 
block of 7 is contained in one of the blocks of a. In this case, we write 7 -< o". Given 
1 ^ i ^ j ^ k, there is one-to-one correspondence between all possible pairs of partitions 
(7, vr) with 7 G Il{k,j) and vr G n(j, £) splitting the set {71, ...,7^} of blocks of 7 into i 
'super-blocks', and all possible pairs of partitions (7,0") with 7 G Il{k,j) and a G Il{k,£) 
such that 7^0". We will use the enumeration ai, ...,ae of blocks of a induced by the 
enumeration of blocks of vr. For instance, in the example considered above, a G 11(8, 2), 
(71 = {1, 2, 3, 4, 5, 6} and 0-2 = {7, 8}. 

We denote by P^,^ = P^^tt, 1 ^ i ^ ^, the set of lengths of blocks of 7 that are 
contained in the block at of a. Then changing the order of summations in (38), we get 

E(-i/-^(^-i)! E E E Ilr-{h,K.) 

i=i aeu(k,i) j=i 7en(fc,j) t=i 

e=i (Ten{k,e) y£n{k) t=i 

7-<cr 

Next, given a G Il{k,i), we fix an arbitrary enumeration of its blocks, put qt = \crt\, 
1 ^ t ^ £, and replace one partition 7 ^ cr by £ partitions 7"*^ G n(gi), ..., 7^ G Il{qe) 
that split the blocks of a into the corresponding sub-blocks. In the example considered 
above, qi = 6, 52 = 2, the partition 7^ consists of two blocks {1, 2, 3} and {4, 5, 6}, and 
the partition 7^ consists of one block. 

At last, replacing the sum y, by the £-tuple sum N^ ... y, , ^^ g^t 

76n(fc) 7^Gn(gi) 7*6n(g,) 

7-<a- 

e e 

7en(fc) i=i i=i Aenfe) 

7-<CT 

Hence, the right-hand side of (33) equals 

i=\ o-6n(fc,£) 

proving the claim. D 
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